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Active diffusiophoresis - swimming through interaction with a self-generated, neutral, solute gradient - is a 
paradigm for autonomous motion at the micrometer scale. We study this propulsion mechanism within a 
linear response theory. Firstly, we consider several aspects relating to the dynamics of the swimming particle. 
We extend established analytical formulae to describe small swimmers, which interact with their environment 
on a finite lengthscale. Solute convection is also taken into account. Modeling of the chemical reaction 
reveals a coupling between the angular distribution of reactivity on the swimmer and the concentration field. 
This effect, which we term "reaction induced concentration distortion", strongly influences the particle speed. 
Building on these insights, we employ irreversible, linear thermodynamics to formulate an energy balance. 
This approach highlights the importance of solute convection for a consistent treatment of the energetics. 
The efficiency of swimming is calculated numerically and approximated analytically. Finally, we define an 
efficiency of transport for swimmers which are moving in random directions. It is shown that this efficiency 
scales as the inverse of the macroscopic distance over which transport is to occur. 
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I. INTRODUCTION 

Active motion, driven by an inhomogeneous chemi- 
cal surface reaction, has recently attracted much scien- 
tific interest. A number of different systems, employing 
this propulsive mechanism on the micro- and nanometer 
scale, have been suggested^"— . They open exiting, pos- 
sible routes for future engineering of nano-devices^rJ^. 
Furthermore, self-propelled small swimmers provide a 
unique chance to study artificial motion at the nanome- 
ter scale, e.g., to investigate anomalous diffusioni^"— or 
non equilibrium thermodynamic o^^'^*^ . Other interesting 
questions relating to mutual interactions of these swim- 
mers and interactions with the environment are also be- 
ginning to be addressed2i"— . 

In most of these experiments, a catalytic decompo- 
sition of hydrogen peroxide H2O2 into oxygen O2 (g) 
and water is employed for propulsion. In general, 
these (electro) chemical processes obey quite nontriv- 
ial kinetics^^. Experiments with bimetallic nanorods, 
combining two different metallic catalysts, have shown 
that the mechanism for electrokinetic decomposition 
of H2O2 involves an electrical current inside the 
particle^i^. Theoretical interpretations of these ex- 
periments in terms of a phoretic propulsion mecha- 
nism have recently emerged22r— . They include driving 
through an interaction of the swimmer with a concen- 
tration gradient (diffusiophoresis^^) as well as driving 
through interaction with a self-generated charge gradi- 
ent (self-electrophoresis). Other possible driving mecha- 
nisms like interfacial tension gradients^ and nano-bubble 
formation^^ have also been investigated. For a different 
kind of experimenti^^i polystyrene particles, half coated 
with platinum as a catalyst, have been employed. No 
electric current inside the swimmer is expected here. 
Michaelis-Menten-like kinetics for the chemical reaction 
were demonstrated and the data can be consistently ex- 
plained by a generic diffusiophoretic/self-electrophoretic 
model^. 



Although the understanding of phoretic driving mecha- 
nisms has reached an appreciable level, there remain open 
questions concerning swimming with progressive minia- 
turization of the swimming particle. One such question 
is, how swimming speed changes if the particle size ap- 
proaches the lengthscale of solute-swimmer interactions. 
Further challenges concern a more detailed, theoretical 
modelling of the chemical surface reaction and the so- 
lute convection. Both factors can contribute to modified 
predictions for the swimming behavior. All these issues 
may also be relevant in terms of optimization for future 
applications (see citations above). 

Finally, the energy balance of the driving mechanism 
and its efficiency have only received limited attention. 
Early experimental estimates^ suggest that the efficiency 
of bimetallic nanorods is very low (on the order of 
10"^). Here the swimming speed was found to be 
about 1 . . . 10 /im/s. Molecular dynamics simulations of 
a swimmer in gaseous environment^ yielded much higher 
efficiencies of the order 10~^ and speeds in the order of 
m/s . The discrepancy emphasizes the importance of 
the hydrodynamic dissipation for the efficiency. This 
point found particular emphasize in our previous the- 
oretical work were we have examined the hydrodynamic 
efficiency for generic surface driven swimmers^. The hy- 
drodynamic efficiency provides an upper bound on the 
overall efficiency while it does not include dissipative ef- 
fects related to building up the chemical gradient driv- 
ing the swimmer. It is a remarkable property of active 
phoresis that the overall efficiency is fully amenable to 
theoretical calculations. This contrasts with other ac- 
tive swimming mechanisms, e.g., beating fiagella^^, where 
some microscopic, mechanical degrees of freedom are un- 
known. Nevertheless, the overall efficiency of a model 
for active, phoretic swimming in a viscous medium has 
not yet been calculated. In the present publication we 
address the above mentioned issues concerning the dy- 
namics and efficiency within a linear response theory em- 
ploying the established theoretical framework of neutral 
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diffusiophoresia^. 

This article is organized as follows. In Sec. we in- 
troduce the model system and the pertaining equations. 
These are non-dimcnsionalized and linearized in Sec. IIIII 
In Sec. IIVI we calculate the speed of the swimmer and 
analyze its dependence on the interaction between solute 
and swimmer. The dependence of the swimming speed 
on the chemical surface reaction is investigated in Sec. 
IVl In Sec. I VII we focus on the energetics of active diffu- 
siophoresis. Here, the efficiency of swimming is analyzed 
in detail and we also consider a measure for the efficiency 
of transport. 

II. MODEL DESCRIPTION 

To explore generic features of active diffusiophorcsis, 
we employ an idealized hydrodynamic model. The swim- 
mer, a spherical particle with radius R, is placed in an 
infinitely large container. We assume axial symmetry 
and use a spherical coordinate system aligned in the 
Bz direction, r is the distance from the particle center 
and 9 is the polar angle with < 9 < tt. The unit 
vectors of the spherical system arc denoted by e^, eg. 
Throughout this publication, wc will designate dimen- 
sional variables with a tilde ( ~ ), while constants and 
non-dimensional quantities carry no tilde. The multi- 
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FIG. 1. Schematic representation of the model for an active 
swimmer. Molecules of type a and b are dissolved in the fluid 
surrounding the swimmer with radius R. The molecules are 
transformed into each other through chemical surface reac- 
tions and the overall concentrations are maintained constant. 
An interaction between the swimmer and the type b molecules 
occurs through the potential whose lengthscale is given by 
L. An inhomogeneous distribution of solutes causes the swim- 
mer to move with the speed U. r is the position vector and 
6 is the polar angle in a spherical coordinate system used 
throughout. 

component fluid surrounding the swimmer is assumed to 
be incompressible and Newtonian. It contains two dilute 



(non-interacting) solutes, whose molecular volumes and 
masses are to be similar to that of the solvent. The con- 
centrations of the solutes are denoted by Ca and Cb- So- 
lutes and solvent are bound to interact with the surface 
of the swimming particle, e.g., through van der Waals, 
electrostatic or steric interactions. In order to simplify 
the situation, we assume that a common surface poten- 
tial mediates the interactions between both, the solutes, 
the solvent and the swimmer. Within our mean-field 
framework, the common interaction between swimmer 
and all fluid constituents merely modifies the pressure 
of the fluid. Only for species b, we assume, on top of 
this, a radially symmetric potential 'I'. Physically, this 
second potential could, e.g., represent a different polar- 
izability of the solutes b. It is assumed to decay over a 
characteristic lengthscale L. A similar approach has been 
taken in a molecular dynamics simulation of the system^. 
In Appendix |A] we provide a simple physical model of 
the three-component solution considered in the present 
publication. Both molecular species, type a and 6, have 
the same diffusion constant D. The concentration fluxes 
J{a,fa} and the corresponding solute conservation laws are 

3a = Ca-V +ja = Ca-V - DVCa, (1) 

3b = £bv + jb = 56V - D (vcb + -^V*^ , (2) 

= y-~3{aM- (3) 

The diffusive fluxes of solutes relative to the center of 
mass are denoted by i{a,b}- The center of mass fluid 
velocity that convects the solutes is denoted by v. We 
assume here that the diffusion of solutes near the surface 
is similar to the bulk diffusion. An extra surface diffusion 
term is not taken into account. Generic boundary condi- 
tions, allowing for an absorption and emission of solutes 
at the swimmer's surface, are 

er3{a.b}{R, 0) = a{a,b} g{0), 
Vc{a^b}(cx),6') = 0. 

Here all the 6'-dependence of the flux at the surface is 
contained in the dimensionless function g{9). We will 
always assume that the emission of solutes happens in 
a spatially asymmetric way. has the dimensions 

of solute flux and quantifies its magnitude. The Stokes 
equation and incompressibility condition read 

r^V^v - Vp = V • o- = CbV* and V • v = 0. (5) 

Here 77 is the viscosity of the fluid, which we assume to 
be spatially constant. The pressure accounting for the 
incompressibility of the fluid is denoted by p and 

CT = 77(Vv + (Vv)'^) -pI = 277E-pI (6) 

is the hydrodynamic stress tensor. In order to make sure 
that the total force on the fluid is, even for a linear con- 
centration gradient (cb ~ fcos9), convergent wc assume 
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in the foUowing 



4' (f ^ oo) < O 



(7) 



The swimming speed is given by U in the laboratory 
frame. In the frame moving with the center of the swim- 
mer, it enters the solution of the Stokes equation through 
the boundary conditions 



v(oo, 
v(i?, 



= 0. 



(8) 



The swimming speed U is determined from the balance 
of forces acting on the swimmer in the ovcrdamped limit 



= F„e, 



(9) 



Since the swimmer is spherical, we can identify the nor- 
mal vector on the surface of the sphere with the unit 
radial vector e^. We are considering a freely moving 
swimmer, therefore the external force -Fm will be set to 
zero. The dependence of a on U allows the calcuation 
of the latter quantity from Eq. (jH]). In a pcrturbativc 
approach for small ^/kT, one can assume that the body 
force V • (T = CbV^ is independent of the fluid velocity v. 
Due to the linearity of the Stokes equation Eq. ([9]) can 
be reformulated employing a reciprocal relational. The 
result for spherical swimmers is 



- 67r?yi?C/ 



3R 

2f ~ 2f3 



1) cos 6 e.rV • (7 



-{^ + ^-l)smeegy ■a]dV. 



(10) 

However, when convection strongly modifies the solute 
distribution c^, Eq. ([TU]) does not permit to calculate U. 

The swimmer model presented above is only valid if 
the assumption of a dilute solution holds. Corrections 
to the diffusion coefficient, resulting from mutual solute- 
solute interactions, are proportional to the volume frac- 
tion and we hence demand c^^.b] <C 1 for solutes with 
radius a. Another source of error is introduced by the 
rotational diffusion of the swimmer. The rotating frame 
results in a corrective term for the solute fluxes Eq. (|ll2p 
of the form — e^r Djot c^gCja.b} where Drot is the rota- 
tional diffusion coefficient of the swimmer. The relative 
magnitude of this correction for the diffusion equations 
is of the order R^D^ot/D. Employing the Einstein 
relations for the diffusion coefficients D = kT/ [Q-K-qa), 
D,ot = kT/{8nr]R^), we have R^D.^t/D = ia/{AR). 
Accordingly, we demand in the following a R. This 
assumption also justifies disregarding hydrodynamic cor- 
rections of the solute-swimmer interaction due to the fi- 
nite size of the solutes. Taken together, our model is re- 
stricted to swimmers which are at least two magnitudes 
larger than the solutes. If the solute radii are in the A 
range, we hence require i? > 30 nm. 



III. THE LINEARIZED, 
EQUATIONS 



NON-DIMENSIONAL 



Convective transport renders the diffusion equations 
(IT][3|) nonlinear since the fluid velocity and the concen- 
tration flelds are both unknown. The fluid velocity v 
is determined by Eq. ^ and thus depends itself on Cfc. 
In order to make the system more amenable to analytic 
calculations we linearize the non-dimensionalized equa- 
tions around equilibrium. We thereby construct a linear 
response theory for active diffusiophoresis. The natural 
length- and energy scales are the radius of the swimmer, 
R, and the thermal energy kT. Accordingly, we rescale 
as r = f/R and = /kT. The intrinsic lengthscale L 
of the potential ^ is denoted in dimensionless form as 



A = L/R. 



(11) 



Particle swimming is driven by an asymmetric concentra- 
tion perturbation of the equilibrium distribution of &-type 
solutes c^{f). Owing to the radial symmetry of ^I', only 
the dipole moment of the concentration perturbation 
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D 



g{d) cos6'd(cos6') 



(12) 



contributes here. We can accordingly define the concen- 
tration scale c, which is relevant for the particle motion, 
as 



c = Rab/D. 



(13) 



In the following we will always assume that the concen- 
tration perturbation c is much smaller than the equi- 
librium concentration scale. We can then define a new 
dimensionless parameter as 



<5 = c/5°''(oo). 



(14) 



Throughout this publication, we employ c for the non- 
dimensionalization of concentrations. The equilibrium 
perturbations of the concentration fields are 



Cb 



Cg - 5gq(oo) 
c 

~cq/ \ — * 

_ Cb - V(oo)e 



(15) 
(16) 



The solute fiuxes j are non-dimensionalizd by Dc/R and 
obey the boundary conditions 

erj{a,6}(l, 0) = ^-j^^ = a{a,b}.9(6')- (17) 

A typical diffusiophorctic speed magnitude^ is given for 
A < 1 by 



U 



kTL"^ c 
T] R' 



(18) 
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Wc define an associated Peclet number by 



Pe = 



UR 



(19) 



Sample calculations, employing measured particle 
speeda^^"^, show that Pe is typically small, on the or- 
der of 10^^. The Peclet number thus constitutes, besides 
5, a second small expansion parameter. Employing the 
non-dimensional hydrodynamic variables v = v/C/ and 
p = pR/{Ur]), the equations for solute conservation ([T][3]) 
and the Stokes equation ([5]) become 



Pe 

V^Cfc + V • (cfcVI') - — vVe"* = PevVcb, 





rCfcV* = 0. 



The fluid boundary conditions are, as before, 

v{oo,0) = -Ue^, 
v(l,0)-O. 



(20) 
(21) 
(22) 



(23) 



The form of Eqns. (|20I21|) suggests a perturbation 
scheme in Pe to cope with the nonlincarities. However, 
difficulties arise if S is of the same magnitude as Pe . Then 
only the right hand sides of Eqns. (|20l2ip are small. Fur- 
thermore, because 5 ~ c and also Pe ^ c, it is in the spirit 
of a linear response theory in c to set 



PevVca ~ PcvVcf, ~ 0. 



(24) 



This linearization approach has the merit that the mo- 
bility of the swimmer, relating U to a^^.b}: can be cal- 
culated straightforwardly. To do this, one only needs to 
determine the force balance Eq. ([9]) for the two cases 
(^{a.b} = and U ~ 0. This procedure was described 
in detail by O'Brien and Whito^. For concrete numer- 
ical and analytical calculations we always rewrite the 
Stokes equation (j22p by employing the stream function 
formalisni22i. 

The approximation Eq. (|24p does not go without any 
caveat. The convection terms dominate the diffusion 
terms in Eqns. (|20l2ip for large r due to different orders 
of radial derivatives. In general, one would therefore fol- 
low Acrivos and Taylor— and split the solution for the 
concentration field into an inner solution, with coordi- 
nate r, where diffusion dominates and an outer solution, 
with rescaled radial coordinate Per, where convection 
dominates. We avoid this procedure by assuming that 
APe <C 1. Then the interaction between solutes and 
swimmers takes place in the diffusion-dominated region 
around the swimmer. 

In contrast to the approach taken here, the standard 
way to make analytical progress in view of the difficulties 
inherent in Eqns. (HJS]) is to assume that the magnitude 
of the surface potential is small (^(1) ^ 1). Then con- 
vection does not occur in the O (\I'(1)^) contribution of 



a regular perturbation scheme. The resulting "Debye- 
Hueckel" -like theory is practical in the case of ionic 
surface interactions^ but may become inappropriate for 
non- ionic interactions when ^"(1) is not small. 



IV. SWIMMING SPEED AND THE ROLE OF THE 
INTERACTION POTENTIAL 

The nature of the potential ^P, comprising a variety of 
possible physical interactions, is important for the diffu- 
siophoretic speed. In order to demonstrate this in this 
section, we leave here all issues concerning the details of 
the chemical reactions aside. This simplification also fa- 
cilitates a comparison with other work^fS. The emission 
rate of the solute of type b is in this section given by 
g{9) = (1 + cos 6*) and ai, = k ^ const. We then have 



(25) 



and c = kR/D. This boundary condition for the surface 
interacting solute does not contain Ca- Species a is thus 
irrelevant for the particle swimming. We hence disregard 
Ca in this section. The concentration field Cb is deter- 
mined by the differential Eq. (^0]) with Eq. ([Ml) . In spite 
of its linearity, there seems to be no analytical solution 
of this equation available for arbitrary ^P. This is even so 
in complete absence of convection when Pe /S = 0. 



A. Analytical approximations 

For short ranged potentials, A ^ 1, one can resort to 
a technique of matched asymptotic expansions to calcu- 
late the concentration field^. We find for the swimming 
speed to lowest order in A, (see Appendix IbI) . 
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(26) 



where j/ = (r — 1)/Aisa dimensionless inner variable 
which is 0(1) near the surface, where ^ changes strongly. 
The unnormalized moments of the equilibrium excess sur- 
face concentration are defined as 



y 



1 dy. 



(27) 



As in the classic electrophoretic Smoluchowski limit^, 
the dimensional speed Uq is found to be proportional 
to the square of the interaction Icngthscale Ki). 
We calculate the first speed correction up to 0(A) for 
active diffusiophoresis by employing methods outlined by 
Anderson, Lowell and Prieve^. The corresponding first 
order Padc approximant of the swimming speed reads 



U^Uo/ 
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Pe M N ~ 



(28) 
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where the following definitions are employed 
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y' (e-*(«') - l) dy' (e-*(^) - l) dy, 



(29) 



TV = -2 



oo 



y' (e-*(^') - l) dy' (e*^'') - l) dy 
(e-*(^) - l) - l) dy. 



JO 

(30) 

Note that the chemical reaction must occur at a finite dis- 
tance from the physical hard-core boundary of the swim- 
mer. Therefore, the potential ^(r — 1) cannot diverge 
in active diffusiophoresis and all the constants defined 
above remain finite. 



B. Numerical results 
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FIG. 2. Swimming speed vs. interaction lengthscale A for an 
exponential repulsion: ^(r) = exp [— (r — 1) /A]. (Dashed 
line) Swimming speed in lowest order approximation A^|l7o 
for A <C 1 from Eq. (|26p . (o) Numerical solution for Pe/(5 = 
0. (•) Fade approximant, Eq. (gS)), for Ve/5 = 0. (□) 
Numerical solution for Fe/(5 = 5. (■) Fade approximant for 
Pe /5 = 5. (A) Numerical solution for Fe /5 = 10. (a) Fade 
approximant for Fe /5 = 10. 



In order to go beyond the A <C 1 limit, we solve 
Eqns. (|20II22I24|) numerically for the boundary condi- 
tions given in Eqns. (|23I25|) . Results are displayed in 
Figs, m [3] and m where the swimming speed is plotted as 
= LP'U / {B?U) in order to demonstrate the physi- 
cal dependence of the speed on the interaction length L. 
As it is the case for other phoretic effects^i^, convec- 
tion causes a non-monotoneous relation between particle 
speed and A (Fig. An increasing magnitude of the 
surface potential ^'(1) reduces the range of validity of 
the lowest order approximation Uq to smaller values of A 
(Fig. El). The Fade approximant of U, Eq. (gSl), then de- 
scribes the case of strong surface interactions much more 
satisfyingly^. Therefore, the Fade approximant presents 
a significant improvement over the lowest order estimate 
f/o, in spite of being useful only for A < 1. 

Fig. m shows the swimming speed for a truncated 
van der Waals-like interaction where the potential de- 
cays ~ l/r® only far away of the swimmer. The result 
demonstrates that multiple lengthscales of the interac- 
tion potential may modify the simple sealing of Uq ^ L^. 
Furthermore, the similarity of the results for Pe/(5 = 
and Pe /(5 = 5 qualifies the simple notion that convection 
of solutes reduces the speed of the particle. Since, for 
Fig. m we have '${r) <C 1 when A > 1, the influence of 
the Peclet number is suppressed almost completely. This 
dependence of the effect of convection on the strength of 
the surface interaction also appears through the constant 
M in the Fade approximant given in Eqns. (|28l29p . 

C. Comparison with passive diffusiophoresis 

It is of interest to compare the swimming speed given 
in Eqns. (|26l28p with the analogous formulae for passive 
diffusiophoresis in an externally imposed concentration 




FIG. 3. Swimming speed vs. interaction potential magnitude 
*(1) for *(r) = *(1) exp [- (r - 1) /A] with Fe /S = 0. (□) 
Numerical solution for A = 0.05. (■) Fade approximant, Eq. 
dlSl), for A = 0.05. (Dashed line) A^jf/o| from Eq. ^ for 
A = 0.05. (o) Numerical solution for A = 0.1. (•) Fade 
approximant for A = 0.1. (Dotted line) A'^|C/o| for A = 0.1. 
Inset: Relative deviation of the lowest order approximation 
{Uo-U)/Uq vs. *(1) with Pe/S = 0. (A) A = 1. (o) 
A = 0.1. (A) A = 0.01. 



gradients^. Here the boundary conditions for the con- 
centration of the solute of type b read 

e.]b(l,^)=0, ^3^^ 
Vc|oo — const. X §2. 

The lowest order result for swimming speed in passive 
diffusiophoresis f/^ is 

C/P = |Vc|oo^^i^i. (32) 
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FIG. 4. Swimming speed vs. A for a van der Waals -like 
attraction: *(r) = -A\^/{r - 1 + \f /{r - 1 + A + 2)^ witli 
A = 1. A is here interpreted as solute radius divided by the 
swimmer radius and the potential is truncated at a distance A 
away from the surface. (Dashed line) A^|l7o| from Eq. (|26l) . 
(o) Numerical solution for Pe /5 = 0. (•) Pade approximant, 
Eq. ((lH), for Pe /S = 0. (□) Numerical solution for Pe /5 = 5. 
(■) Pade approximant for Pe /6 = 5. 



The formulae ((26)) and ([32)) for active and passive swim- 
ming are, apart from a replacement of k/3D by |Vc|oo, 
the same. These different prefactors result from the 
dipole moments of the concentration fields around the 
swimmer. The concentration far away from the parti- 
cle is given by Cf, « Rk/D [l/r + cos6'/(2r^)] in our case 

and Cb w i?| Vc|oo [r + 1/(2^'^)] cos 9 for an externally im- 
posed concentration gradient. At r = 1, the self gener- 
ated dipole is therefore only ^ kRD/2 while the imposed 
concentration dipole is 3i?|Vc|oo/2. This leads to the 
differing factor of 1 /3, appearing in the swimming speeds. 

The general agreement of both lowest order formulae 
for active and passive diffusiophoresis can be rationalized 
by noting that both concentrations are, in close proxim- 
ity to the swimmer's surface, in radial equilibrium. Ac- 
cordingly, the lowest order solute flux also vanishes for 
active diffusiophoresis since the diffusive exchange of so- 
lute near the active swimmer is ^ dcb/d{Xy) = 0(1/A) 
while the emission/ absorption of solute, determined by 
g{e), is only of 0(1). 

For passive diffusiophoresis^, the analogous formula 
to the first order swimming speed Eq. ()28p is, in our 
notation. 
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K2 Pe M 
2KI ^ TY 



(33) 



The main difference between (|28)) and (|33)) is the ap- 
pearance of the term XN/Ki for active diffusiophoresis. 
The constant N, defined in Eq. (jSO)) . contains an in- 
tegral over exp(4'(r)). It increases the relative impor- 
tance of the correction for strong repulsive surface 
interactions. For monotonic potentials, iV is a positive 
quantity. Therefore, the correction for active swimmers 
increases the swimming speed for purely repulsive inter- 
actions where we have Ki < 0. Convective corrections 



^ Pe /S are, both for active and passive swimming, rel- 
evant if the equilibrium concentration scale c^'^ is much 
larger than the concentration disturbance c. It is an in- 
teresting side note to the convective correction that the 
constant M is related to the hydrodynamic dissipation 
in the boundary layer^ Whyd « M An {kT cf / {3r]). 



V. REACTION INDUCED CONCENTRATION 
DISTORTION 

In general, the kinetics of the chemical transformations 
at the swimmer's surface depend on the local concentra- 
tions and correlation functions of educts and products. 
An inhomogeneous reactivity profile on the swimmer, re- 
alized, e.g., through partial coating with a catalyst, cou- 
ples to the angular solute distributions. Therefore, the 
resulting solute flux at the swimmer will usually not have 
the same angular dependence as the reactivity profile. 
We term this effect "reaction induced concentration dis- 
tortion" and it bears a certain analogy to polarization 
phenomena, inasmuch as it has a possibly nonlinear ef- 
fect on the swimming speed. 



A. Boundary conditions for the concentration field 

For the considered case of dilute solutions one might 
anticipate that the fraction of unoccupied catalytic sites 
on the swimmer is always very small. Then the overall 
reaction rates depend linearly on the concentrations of 
the reactants. Concentration fields resulting from non- 
linear reactions^ are beyond the scope of this article. 
We therefore assume, for simplicity, a reversible, first or- 
der reaction of the form a -O- 6. For the remainder of 
this article we will employ as boundary conditions for 
the concentration fields at r = 1 

er]b{l,0) = (1+COS0) [khCa{l,0)~hM^,0)\ , (34) 

e.ja(l,e) = -6rjh(l,e), (35) 

where fcab and fcba are rate constants which we non- 
dimensionalize with D/R. The concentration depen- 
dence of the flux at the swimmer surface is also termed 
a radiation boundary condition. Similar mathematical 
problems have occurred in the calculation of mean first 
passage times for the combined encounter and reaction 
of asymmetric molecules^i^. 

In equilibrium, where the fluxes vanish, the radially 
symmetric distributions of solutes obey c^{r) = d^loo), 



Hoo)e-*('^) and 
(00) 
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(00) 



k 



(36) 
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A flnite reaction rate at the swimmer's surface is ul- 
timately driven by chemical potential differences far 
away, at the boundary of the system. Expressions 



7 



for the chemical potentials are given in Appendix |^ 
Within the linear response regime we expect that 

(^C{a,6}(oo)/c^^ ^j(oo) — 1^ <^ 1. Therefore, the chemi- 
cal potential difference at f ^ oo becomes 



ilgjoo) - /2b(oo) _ Cg (oo) ^ Cb (oo) 



kT 



1 



The concentration scale Eq. (|T3|) can now be defined as 

(38) 



R 



and from Eq. (g]) with Eq. ([M]) we have g{6) = 
erjbll,^') X R/{Dc). The parameter (5 and the velocity 
scale U arc accordingly given by 



1/ 



77 R 



(39) 
(40) 



In order to calculate the concentration perturbations 
C{a.b} around the swimmer we proceed by expanding 
them in Legendre polynomials P„(cos0) as 



An 



A^. Finally, Eq. (|42p yields a recursion 



equation for the constants 



2{n + l + k+) 2nk+ ^ 



2n+l 
2{n+l)k+ 
(2n+ 1) (2n + 3) 

Here we defined 



(271 + 1) (2n - 1) 



(44) 



fc+ = kah + khi,e = fcab ( 1 



'(oo) 



'(oo) 



(45) 



Ea. (|44)) can be written in matrix form as Bj = Mjn A^. 
The off-diagonal elements of {Mj„} decay like ~ l/n 
and we can invert a finite matrix to determine a 

numerical approximation of the {^^'}. In the following, 
plots of analytical results involving are created by 
employing Eq. (|44| and rimax = 70. Since \E'(r) is ra- 
dially symmetric, the swimming speed depends only on 
the dipole moment of the concentration field. Accord- 
ingly, an expansion of Eq. (|T0)) for small A yields for the 
lowest order the free swimming speed 



C/n = 2Al 



(46) 



Ca ^Ca (oo) + ^ c2{r)Pn {cOS 6) , 



(41) 



Cb =Cb 



(oo)e-*M+^cnr)P„(cos0), 



n=0 



which includes the boundary conditions at r — >■ oo. The 
boundary conditions at the surface of the swimmer, Eqns. 
P4l35p , couple different coefficients of the expansion Eq. 
(pi and one has 

[drc'^ + dj^dr"^] \r=i = / P„ (cos 6*) g(6') sin6'd6'. 



2n+ 1 



B. Analytical approximations 



(42) 
(43) 



We consider a swimmer with short interaction length 
A <C 1. Employing the methods presented in Appendix 
IB] we calculate the concentration fields to leading order 
in A. Far away of the swimmer, where \E'(r) 0, Eqns. 
(|20|24p can be replaced by Laplace's equations and we 
have c7^f,i(r) = A? in Eq. (gT]). For Ca, this 



expansion is valid throughout the whole system. For the 
6-type solute we have to leading order near the surface of 
the swimmer c^{y) « a^e~*'^^. Matching these solutions 
and employing the boundary condition Eq. (j43|) leads to 



For small bare rates, and therefore small one can 
expand the inner concentration field of the 6-type solute 



Cb ==Cfc(oo)e-*W+e-*Wx 

[^^ + Al cos 9 + AlP2{cose) + AlPsicosl 



+ 0{kl) 
(47) 



with following approximations for the constants 



^b 



^b 



42 + 2k+ 
42 4- 51k+ 
2k+ 



18 + 33fcj 



41 

^b 



43 
^b 



135 - 2Qk+ 
270T353fc^ 

K 
60' 



(48) 
(49) 



The magnitude of the dipole moment is determined by 
A\ Ki (1/2 — 3fc+/4). Reaction induced concentration 
distortion, emerging here through the corrections in or- 
ders of fc+ , reduces the dipole moment and thus slows the 
particle swimming down. 



C. Numerical results 

To complement the analytical approximations, we cal- 
culate the swimming speed numerically. Eqns. p4l35p 
with fixed concentrations far away of the swimmer are 
employed for numerical solution of Eqns. (|20ll24p . The 
boundary conditions result in an infinite system of equa- 
tions where the solute concentrations far away of the 
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swimmer C{a^f,}(oo) determine the reaction speed. The 
system can be truncated above a certain order nmax of 
the Legendre polynomials. In choosing here rimax = 4 we 
made sure that the error in the calculated concentrations 
is negligible. Fig. [5] contains an exemplary plot of how 
the reaction induced concentration distortion influences 
the swimming speed U. The chosen velocity scale U (Eq. 
HO)) contains the linear dependence on C{a,b}{oo). With 




FIG. 5. Particle speed [/ as a function of fcab and fcba for 
vE'(r) = exp [- (r - 1) /A], Pc/5 = 1 and A = O.f. 

the employed first order reactions, the concentration dis- 
tortion depends nonlinearly on the bare rates fcab and fcba 
but linearly on the concentration scales C{a,b}(oo)- Thus, 
in an experiment in the linear regime, with fixed fcab and 
fcba, the reaction induced concentration distortion might 
be accounted for by a constant prefactor, modifying the 
swimming speed. 

lUl 

0.30t- 

0.25 



0.20- 




FIG. 6. Rate constant dependence of the swimming speed 
U for ^(r) = exp [- {r - 1) /A]. See Eq. ^ for the def- 
inition of k+. (Dashed line) Reaction induced concentra- 
tion distortion neglected as explained in the text, leading 
to ;7 = Ki/3 with Ki ~ -0.89. (Dotted line) Approxima- 
tion Uo (Eq. (|46p ) for A ^ f which takes the concentration 
distortion into account. Full symbols are numerical results 
for fcba = and open symbols are results for fcab = 0. (♦) 
A = 0.005, Pe/5 = 0. (■,□) A = O.f, Pe/S = 0. (.,0) 
A = O.f, Pe /5 = f 0. (A, A) A = f , Pe /5 = 0. 



In Fig. |6]we plot numerical results for the swimming 
speed as a function of defined in Eq. dJS]). Ne- 
glecting the reaction induced concentration distortion, 
the naive boundary condition for the concentration would 
be erjf,|r=i = (1 + cos9). The resulting swimming speed 
is independent of k^^. It agrees with the analytical ap- 
proximation in Eq. (|46p in the limit — >■ 0. However, 
Fig. [H] shows that neglecting reaction induced concentra- 
tion distortion in this way leads to significant errors in 
the speed estimate for finite reaction rate constants. The 
analytical approximation, Eq. (j46p . is found to be useful 
for A < 0.01. 

Fig. [7]shows the dependence of swimming speed on the 
strength of the interaction potential ^'(1). The symmetry 
between the effects of changing fcab exp (— 'I'(l)) and fcba, 
apparent in the parameter fc+ in Uq when A <C 1, is lost 
for finite A. Increasing the strength of the interaction 
potential ^'(1) makes this asymmetry more pronounced. 

U 




FIG. 7. Swimming speed as a function of the potential 
strength ^'(1). Numerical results are computed with 'l'(r) = 
*(f)exp [- (r - f) /A] and Pe/5 = 0. Full symbols indi- 
cate (fcab ~ f, fcba exp (^^(f)) — 0). Open symbols indicate 
(fcab = 0, fcba exp (*(f)) = f). (^,0) A = 0.005. (■,□) 
A = O.f. (A, A) A = f. 



VI. ENERGETICS 

A. Energy balance 

We consider an isothermal situation (T = const.) 
where the system, i.e., the swimmer and the multi- 
component fiuid is in a steady state. The steady state 
implies, that the molecules which are modified by a chem- 
ical reaction need to be replenished from outside the sys- 
tem. To do this in a real experiment, one would need to 
connect some sort of external apparatus to the system. 
We idealize the apparatus by a reversible process con- 
suming the power Pin. The total of external apparatus 
and system does not exchange matter with the exter- 
nal world. Therefore Pin is balanced by the overall heat 
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outflow from the system and apparatus. On employing 
classical, linear nonequilibrium thermodynamics wc find 
the power input (see Appendix [C]) 



Pin = kTcDR J (Ohyd + Odiff + e,eact) dV. 



(50) 



The power is consumed by the following three entropy 
production rates per volume 6(...) 



©hyd = -Pc^ 



Cb 



(51) 



©diff = -jfcV* - (52) 

i 

©react = <S (r - 1) gi9) (Aia - * - Mb) • (53) 

The integral over ©hyd is the overall hydrodynamic work 
because no external force is applied to the swimmer and 
because V • v = 0. Its proportionality to Pe emphasizes 
an interpretation as convection of 6-type solute within 
the potential 4". Furthermore, convection of solutes also 
modifies the reaction rate through the mass action law 
underlying the definition of g{0). Hence, the common ne- 
glect of convection in dynamical problems^ would make 
an overall energy balance within the chosen framework 
impossible. In the linear response regime the power con- 
sumption up to 0{6'^) can be rewritten for the employed 
model as (sec Appendix [C|) 



Pi„ = kTDcR / JhA/iooBrdAr 



(54) 



Note that A/ioo is the chemical potential difference at 
the outer boundary of the system, not the local chemical 
potential difference at the the reaction site. Eq. (|54|) is 
nothing but the free energy exchange between the system 
and the apparatus. Due to the similarity of masses and 
sizes of the fluid constituents, it does not make a differ- 
ence whether we consider a fixed pressure at the outer 
boundaries of the system or a fixed system size. 



B. EfFiciency of swimming 

Our freely moving swimmer does not have an external 
power output which could be employed to calculate its 
efficiency. Still, one might ask how efficient this swimmer 
can transport itself. A natural way to do this, is to com- 
pare the energy dissipation of active swimming with the 
energy dissipation taking place when dragging the same 
particle. We accordingly define a swimming efficiency as 



e = ; = PcA 



Pin 



Pe 



«bae 



-*(i);^2 



(55) 



where we have used Eq. (|39p in the second line. The 
numerator of Eq. (j55p is the hydrodynamic dissipa- 
tion of a passively dragged sphere and the denomina- 
tor is the power consumption of our external apparatus. 



providing the energy for active swimming. This defini- 
tion is a natural extension to Lighthill's formula for hy- 
drodynamic efficiency^. The power consumption Pin is 
bounded from below by the hydrodynamic dissipation 
2rj j 'El : Vv dV". This intuitive result follows formally in 
Eq. ([50| from the linear force-flux relationships. Since^ 
GnriR < 2j] J E : Vv dV^, wc always have e < 1 as long 
as no approximation is used to evaluate Eq. ((55)) . 



C. Analytical approximation for the efficiency of swimming 

In the limit A <C 1 one can employ the lowest order 
speed Uq, Eq. (|46)) . for the calculation of the efficiency 
from Eq. ([55]) . The result, including the effect of reaction 
induced concentration distortion, is 



e Pc 



47r.40AMo 



Ds L'^c kT An (AlKi) 
~D' A/2oo 



(56) 



where we have employed the translational diffusion con- 
stant of the swimmer Z?s = kT / Q'K'qR. The first three 
dimensionless groups in the second line of Eq. ([55)) de- 
termine the magnitude of the efficiency. As discussed in 
our previous work^, (nano-) swimmers with interaction 
lengths comparable to their size can have a higher effi- 
ciency than swimmers with L R. This is is evident 
from the factors Dg/D and cL'^/R = cL^ x A in Eq. 
(|56)) . We evaluate Eq. ([56)) by employing the parameters 
Al,Al calculated from Eq. ([33)) with yi^ax = 70. The 
results agrees with numerical solutions for Eq. ([55)) (see 
below) when A < 0.01. However, ignoring convection in 
the denominator of Eq. ([55)) due to the small A limit 
implies neglecting hydrodynamic dissipation. Therefore, 
the asymptotic efficiency in Eq. ()56|) is not strictly 
bounded by unity. 



D. Numerical results for the efficiency of swimming 

For a numerical evaluation of Eq. ()55p we truncate 
the expansion in Legendre polynomials at rimax = 4 as 
in Sec. |V] Due to the linear response nature of our the- 
ory, numerator and denominator of Eq. (|55p are both 
quadratic in A/ioo- Therefore, the equilibrium perturba- 
tion driving the motion of the swimmer does not appear 
in the efficiency. However, the results support the no- 
tion that the swimming efficiency increases away from 
the quasi-equilibrium limit. 

As seen in Figs. [S) and [S) the swimming efficiency is 
proportional to A^ Pe /5 ^ for A 1 . This scaling 
is also evident from the prefactor in Eq. ([55)) . With a 
fixed Peclet number, Pe/5 > 0, the swimming efficiency 
decreases for A > 1 (Fig. [5)). Comparison with Fig. 
[5] shows that the reduction in efficiency is due to the 
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reduction of swimming speed in this regime. When, for 
A > 1, the swimming speed does not increase ^ , 
fixing the Peclet number in Eq. (|55|) introduces a scaling 
of e with a negative power of A. As an alternative, one 
could remove the dependence of Pe on the interaction 
lengthscale L by setting Pe/((5A^) = const for Fig. HI 
This way of plotting the data would render the decrease 
of e for A > 1 less pronounced. 

According to Eqns. (|19I39I40[) . fixing Pe/i5 and A im- 
plicitly sets an absolute equilibrium concentration scale. 
Therefore, Fig. IHlalso suggests that the efficiency of diffu- 
siophoretic swimming increases with the absolute concen- 
tration scale for A < 1. Only the asymptotic analytical 
result for A = 0.01 has been plotted in Fig. [9] because the 
curve for A = 0.1 already showed significant deviations 
from the numerical data. 

Fig. [10] shows the dependence of the swimming effi- 
ciency on the strength of the interaction potential ^'(1). 
\E'(1) influences both, the reaction rate and the swim- 
ming speed and therefore it has a nonlinear effect on e. 
Reaction induced concentration distortion plays an im- 
portant role for the shape of the curve, in particular for 
^'(1) < when A;+ can become much larger than unity. 
For \E'(1) the swimming efficiency vanishes. 

Figure [TT] shows the swimming efficiency for fixed 
equilibrium constant (?^{po) /cj^ipo) and varying reac- 
tion rate /cba. For simplicity, we consider here only 
A = 0.01 <C 1. Due to the truncation of the full 
numerical solution at Legendre polynomials of the or- 
der Timax = 4 the error in the numerical data becomes 
large beyond the plotted range. The reaction induced 
concentration distortion again explains major features 
of the plotted curves. For fcba Gxp(— ^'(1)) ^ 1 and 
fcbaexp(— ^'(1)) < S^^{co) /cj^ {oo) the reaction induced 
concentration distortion is negligible. Then e rises lin- 
early with fcba and the curves fall onto each other. For 
larger equilibrium constants with c°^{oo) /<fj^{oo) ^ 10 
the location of the maximum of swimming efficiency be- 
comes becomes independent of the equilibrium constant. 
On employing Eq. (|56p we here find the maximum of e 
at fcba ^ 1.13cxp*(l). 



E. EfFiciency of transport 

The definition of swimming efficiency in Eq. (j55p is 
based on the comparison of energy dissipation per unit 
time. In certain practical applications of diffusiophoresis 
one might prefer other efficiency measures. It is, e.g., 
interesting to have a measure for the energetic cost of 
using active swimmers for transport between two loca- 
tions. Then, it may be more appropriate to compare 
energy dissipation per transport distance. For concrete- 
ness, we think here of a slab geometry which the swim- 
mers are to cross. Thereby, they move from one slab 
at X = to the other slab sX x = X > Q. Since 
the swimmers are freely suspended, they will not swim 
straightly but Brownian, translational and rotational. 




10"'4=:i; . . . I ....I A 

0.01 0.1 1.0 10 

FIG. 8. The swimming efficiency e as a function of A with con- 
stant Pe /5 for ^'(j/) = exp [—(r — 1)/A] and fcab = ^ba = 1. 
(Dashed lines) Asymptotic results calculated from Eq. (|56|l . 
(A) Pe/(5 = 10. (0) Pe/5 = 1. (□) Pe /5 = 0.1. (o) 
Pe /5 = 0.01. 
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FIG. 9. The swimming efficiency e as a function of Pe /5 
with constant A for ^(y) = exp[— (r — 1)/A], fcab = fcba = 1. 
(Dashed line) Asymptotic result for A = 0.01 calculated from 
Eq. (dSl). (■) A = 10. (A) A = 1. (o) A = 0.1. (A) 
A = 0.01. 



diffusion comes into play. The translational diffusion has 
a constant _Ds = kT / {6Tn]R) and the rotational diffusion 
happens on a timescale Trot = SnrjR^ /kT. For long times 
t ^ Trot J one can employ an effective translational diffu- 
sion constanti^ii^i^i^ where the effects of random diffu- 
sion and active, translational motion are incorporated as 
Doff = Ds + f/^Trot/6. This approach neglects possible 
modifications of Ds and Trot due to the active processes. 
In our setup we expect the transport distance X to be 
much larger than a characteristic length of the random 
walk X ^ Dcs/\U\. Employing reflecting boundary con- 
ditions at a; = 0, the mean first passage time tmf to reach 
X = X is simply given by imf = X^/ (2-Dcff)- The en- 
ergy consumed by the swimmer during the time tmi can 
be estimated by Pin^mf- For an energetic comparison, 
one might consider dragging a passive particle directly 
across X. The corresponding dissipated work would be 
6Tn]R\U\ X . Hence, we define a new efficiency of trans- 
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FIG. 10. Swimming efficiency e vs. the strength of the inter- 
action potential *]>(1), where = >1/(1) exp [-(r - 1)/A] , 
fcab = fcba = 1 and Pe /5 = 0.1. (Dashed line) Asymptotic 
result for A = 0.01 calculated from Eg. ((5^ . (A) A = 1. (o) 
A = 0.1. (A) A = 0.01. 
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FIG. 11. Dependence of the swimming efficiency e on the 
reaction rate constant fcba with fixed equilibrium constant 
= c:''(oo)/c^'*(oo) and Pe /5 = 0.1, A = 0.01, ^{y) = 
exp [— (r — 1)/A]. (Dashed lines) Asymptotic results calcu- 
^ ' 10-^ (o) 



lated from Eq. ((56}. (A) = 10~". (■) 
= 1. (□) = 10^ (A) A'a = 10". 



port as 



etr 



6Tn-iR\U\X _ ^ X 



\U\t„ 



(57) 



which we relate on the right hand side to the efficiency of 
swimming e through Eq. ([55]). On employing the result 
for the mean first passage time i,„f given in the text above 
we find 



et 



ransp 




(58) 



This dependence ~ 1/X of the relative energetic cost 
of direct transport compared to enhanced diffusion is a 
quite generic result. It emphasizes the necessity to im- 
pose a directionality on the active motion if the transport 
distance is to be of macroscopic size. 



VII. DISCUSSION AND CONCLUSION 

In the present publication we take a close look at the 
dynamics and efficiency of diffusiophoretic swimming in 
the linear response regime. In Sec. IIVI we analyze the 
role of the surface interaction potential with regard to 
the dynamics of the swimmer. If the lengthscale of the 
interaction potential is not negligible compared to the 
swimmer radius (A > 0.01), the details of the interaction 
become important. Analytical corrections in powers of A 
to the lowest order swimming speed include various mo- 
ments of the solute concentration and also the effect of 
solute convection. These corrections occur for active dif- 
fusiophoresis in a qualitatively similar way as for passive 
diffusiophoresis^. A particular feature of active diffu- 
siophoresis, however, is that the emission of solutes can 
cause an increase of speed for repulsive surface interac- 
tions. The enhanced speed reflects the effect of "produc- 
ing" solutes directily inside the region from which they 
are repelled. 

In Sec. |V]we employ a linear reaction scheme to model 
the active conversion of solutes at the swimmer's surface. 
Concentration dependent boundary conditions result in 
a strong modification of the concentration fields as com- 
pared to fixing the fiux of solutes at the boundary. This 
reaction induced concentration distortion also occurs in 
the limit A ^ 1 and has a pronounced effect on the swim- 
ming speed. In the linear response regime it causes a 
reduction of the swimming speed, independent of the ab- 
solute concentration levels. If the reaction rate becomes 
nonlinear in the concentrations, a measurable, concen- 
tration dependent, speed modification is to be expected. 
The whole effect may vanish only in the strongly reaction 
limited regime of very high substrate concentrations. 

In Sec. IVII we present results concerning the efficiency 
of diffusiophoretic swimming. The active motion con- 
sidered in the present publication is energetically quite 
distinct from classical phoresis in external gradients^ . 
For the latter systems, the absolute level of concentration 
changes when the particle moves in an externally applied 
gradient. In our case the absolute concentrations remain, 
on average, the same; which implies, that our system 
can be held in a true steady state. The studied model 
includes three modes of entropy production namely hy- 
drodynamic dissipation, dissipation in the chemical re- 
action and entropy production through diffusion of so- 
lutes. While the former two can be recognized imme- 
diately as contributions to the power input, the latter 
becomes only meaningful when considering an external 
apparatus, which converts free energy exchange at the 
outer boundaries of the system into work. Identification 
of the power consumption with the free energy exchange 
at the outer boundaries of the system is a consequence 
of the spatial modelling of the process. It is a natu- 
ral extension to the commonly employed formalism for 
biomolecular motors^, where the local free energy ex- 
change rate is taken as work input. It is pivotal for the 
energy balance to recognize of the importance of convec- 
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tion. For illustration, we think of a situation where the 
reaction rates would be fully independent of the fluid mo- 
tion. This would mean that the net exchange of solutes 
with the external apparatus would be independent of the 
motion of the particle. Thus, according to Eg.lfM]). the 
power input would be independent of the hydrodynamic 
dissipation, which violates energy conservation. Within 
this model, the only way how the particle speed may feed- 
back on the chemical reaction is via a convective modi- 
fication of the concentrations. This link leads, together 
with the mass-action law, to a proper energy balance. 

Relying on the robustness of the scaling of the efB- 
ciency e ~ Ds/D x L'^cjR (Eq. (|56l)). one can employ 
experimental values for catalytic swimmers to estimate 
their efficiency. For swimmers with size in the 1 [im 
range, we have Ds/D ~ 10"'*. Given a concentration 
scaled of c ~ IQ^l/^m? and an interaction length of 
L ~ 1 nm we have L'^c/R = 10"'"'. Together we find 
e ~ 10"^. This number corresponds to the experimental 
estimate for the swimmers of Paxton et. al^. The agree- 
ment may, however, be incidental since these swimmers 
are operating mainly with a self-electrophoretic mecha- 
nism instead of the difFusiophoretic mechanism studied 
here. Naturally, if the concentration perturbation c is 
caused by a large chemical potential difference, this may 
also affect the efficiency as seen from Eq. (|56p . 

In conclusion, we expect the main relevance of this 
work to lie in the scaling predictions and in the iden- 
tification of the main trends in the numerical data. A 
more detailed modeling of specific phoretic mechanisms 
is challenging, but desirable. One interesting question is 
to replace the conservative potential by more realistic, 
angular dependent surface interactions. In reality, the 
solute-swimmer interactions may even transfer energy be- 
tween the fluid and microscopic degrees of freedom on the 
swimmer's surface. Thus, they can acquire a "dissipative 
nature" , possibly resulting in an effective modification of 
the viscosity around the swimmer. Also, lateral diffusion 
of solutes adsorbed to the surface of the swimmer could 
be taken into account here. Issues like hydrodynamic 
surface slip*^ may further complicate the situation. For 
a full description of the behavior of catalytically driven 
swimmers one must also go beyond the linear response 
theory. This refinement of the model will particularly 
effect the reaction induced concentration distortion. Fi- 
nally, the predictions for the scaling of the swimming 
speed and efficiency offer the possibility of being exper- 
imentally tested, which could add interesting new facets 
to the picture of diffusiophoresis presented here. 



The volume variable is written as V . We define a local, 
position dependent free energy per unit volume / = F /V 

as 

/(Ca, Cb, Cs) = 0*0 + CfcUfe*! + CfaCfe + Cata + (0 ^ 1) + 



kT 



Ca log 



h log 



V 4> 



Cs log 



(Al) 



where cj) = CaVa + c^Vb + CsVs- Here the surface interac- 
tions between the fluid constituents and the swimmer are 
incorporated through the potential energies per volume 
^1^0 and ea,eb are the internal energies of the solutes of 
type a and b respectively. The local, hydrostatic pressure 
is calculated as 



P = -F: = -J + -^Ca + -—Cs + TT^C^ = 



dV 



dCa 

W - 1] 



dcs 



dcs 



(A2) 



In the incompressible limit one has cj) = 
(caVa ChVh CsVs) « 1 and the pressure can be 
expanded as p w ;\^[(/) — 1]. The generalized chemical 

potentials /i,^ = /i^ + vi {^q -I- (5i,b*I'i^ = df /dci become, 

with Call = {Ca + Cb + Cs), 



jl's ~kT [log (CsVs) - CallVs] + 

Vs (*0 +p)+ kT {Vsdall - 1) ( 
p.'^ KikT [log {CaVa) - CallVa] + £« + 

Va {^0 +Pj +kT [VaCall - 1) 
KikT [log (cfeUfc) - CallVb] + £6 + 

Vb f *o + ^i+p] +kT {vbdaii - 1) (0 - 1) 



1), 
-1), 



(A3) 
(A4) 
(A5) 



We now apply the highly simplifying assumption of sym- 
metric solute constituents 



Va ^Vb~Vs, 

ma^mb^ rus 



(A6) 
(A7) 



through which particle and mass fluxes become equiva- 
lent. Also, we assume that the volume fraction of solvent 
is much larger than that of the solute 



CaVa + CbVb < 1- 



(A8) 



Appendix A: Simple model for the dilute solution 

In this Appendix we provide a simple, but detailed 
derivation of the model fluid employed in the main part 
of the text. We denote the molecular masses of the so- 
lutes and solvent by ma,mb and mg. The molecular vol- 
umes arc constants, denoted by Va,Vb and Vs respectively. 



The linear phcnomcnological equations, linking a diffu- 
sive flux j to gradients in chemical potential, read 



jb — LbbV 
Ja — Laa^ 



K 




rub 


rrisy 






ma 


rris 



LbaV 
-Lab'^ 







yma 


ms^ 


(K 




\mb 


rUs 



(A9) 
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The flux of the solvent can then be caleulated from the 
condition 



^ mji = 0. 



(AlO) 



■i— a,6,s 



The coefficients Lmn in Eq. (|A9|) relate gradients in 
chemical potential to dissipation. For CqWq + ChVf, <C 1, 
the diagonal coefficients Laa and are proportional 
to the product of a solute volume fraction and the sol- 
vent volume fraction. The solvent volume fraction is 
close to unity and we set L{,b = —rrihCbD/kT and Laa = 
—maCaD/kT. The cross-coefficents, on the other hand, 
are proportional to the product of the two solute volume 
fractions. They obey the Onsager relation Lab = Lba- 
To lowest order in solute volume fractions we therefore 
have Lab = Lba ~ when Eq. (|A8[) holds. On employing 
Eqns. (|A6IIA9p we finally obtain 



ja = -DVca, 

h = -D (\/Cb + ■g;V* ) , 



(All) 
(A12) 



where the surface potential acting on each molecule of 
type b is defined as ^ = Vb^i- Due to the assump- 
tion of similar specific volumes of the molecules in Eqns. 
(|A6IA7p the pressure becomes a purely hydrodynamic 
quantity. We have in steady state the common incom- 
pressibility condition V • v = since 

=V • [{maCa + rribCb + m^Cs) v] « 

TOS_ r,^ ^ . s -1 ™S„ ~ (A13) 

V • [{CaVa + CbVb + CgVs ) vj W V • V. 

Vs Vs 



Appendix B: DifFusiophoretic swimming for A <C 1 

In this Appendix the speed Uq of active swimmers is 
calculated analoguous to work by Anderson and Prieve^ 
for passive swimmers. The diffusion equation (|20l24p be- 
comes in spherical coordinates 



a2c(r, e) + (^^'{r) + drc{r, 9) 



+ ( + -*'(r)) c{r,e) + —^de{smedgc{r,e)) 

r I sm ti 



Pe 



(Bl) 



Due to the radial symmetry of ^(r), only the dipole mo- 
ment of the concentration field can move the particle. 
In absence of a coupling between different spherical har- 
monics it is therefore sufficient to consider the field con- 
tributions which are ~ cosO. We split the concentration 
field of type h solute into an inner field C"(y)cos6' and 
an outer field C°(r) cos 6*. C^{y) is written in terms of 
the inner variable y = {r ^ \) /A. The outer field lies in 



the region with r ^ A -I- 1 where the effect of the surface 
potential ^ is negligible. Eq. (|Bip thus becomes here in 
the outer region 

dlC°{r) + \c°{r)-^C°{r)^Q. (B2) 

For active diffusiophoresis, the concentration perturba- 
tion must vanish far away of the swimmer, Therefore, 
the solution of Eq. (|B2|) is C°{r) = Ajr^ , where the co- 
efficient A needs to be determined by matching with the 
far field behavior of C*. The smallness of A suggests an 
expansion of the inner and outer fields as 



= + C\{y)\ + C\{y)\^ + 0{\^), (B3) 

C° - ^ (Ao + AiA + A2\^ + ...)= Ao+ 

(^1 - 2Aoy) A + (^2 - 2Aiy + 3^oy') + ©(A^). 

(B4) 



In the inner region we substitute r by Ay + 1 in Eq. (jBip 
and expand for small A. The resulting equations for the 
two lowest coefficients Cq and C\ are 

5,(5,CS(y)+CS(2/H*)=0, (B5) 

(B6) 

The boundary condition Eqns. PTI [^5]) yield for the 
coefficients of C* (y) 



-A" 



(a,C;(y)4-C(yHvI;) j^^^ = 



(B7) 

when n^\. 
(B8) 



Employing the above equations, the concentration is 



(e*(^')-l) dy'-y 



(B9) 



The O (A) term diverges for y ^ oo. In order to make 
this divergence explicit, we have removed it from the in- 
tegral in Eq. (jBOp by subtracting 1. The unknown con- 
stants in the inner and outer solution are determined by 
matching them asymptotically^ through C° (y ^ 0) = 
(y ~^ cxd). The resulting conditions, which must be 
valid for all y, are 

^0 = flo, (BIO) 
Ai - 2A„y = ai - H fe**^') - l) dy' - y. (BU) 



This yields the lowest order coefficients aq = Ao = 1/2 
and the innermost concentration field is thus given by 
C wexp(-*)/2. 
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In order to calculate the fluid flow near the surface of 
the swimmer, we employ the Stokes equation (j22[) and 
assume that the body force vanishes in the outer region. 
Applying the curl to Eq. ([^^ and defining a stream 
function S{r) via v = V x {sin 9 S{r)/re^) the Stokes 
equation becomes 



Ad^S{r) 8drS{r) 8S{r) _ C 



4 =-^drnr)- 



(B12) 



In the outer region, where the body force vanishes, we 
have the stream function 5° = X/r + Y r + Z r'^. The 
constants X,Y and Z are expressed as a power series of 
A. The fluid velocity in the outer region thus is 



sin 6 r 
cos 6 



2Zo -Yo + Xo + O (A) , (B13) 
2 ~2Zo + 2Yo + 2Xo + (A) . (B14) 



In the inner region, we expand Eq. (|B12[) for small A 
and insert C"(y) for the concentration field. The leading 
order differential equation for the stream function near 
the surface of the swimmer S*' reads 



dp\y) = -—dy^iy). (B15) 



This yields 



S'iy) ==fco + /o y + Too 2/^ + no - ^h{y) ^^^^^ 

A (fci +liy + mi y^ + m y^) + O (A^) ; 
py rv' , s 

y ^^ (e-*(^ ) - 1) d/vv- (B17) 

The fluid flow in the inner region becomes 

~ -T- (^0 + 2TOoy + 3?io2/^) + \ 'dyh{y)+ 
smt' A ^ '2 

{-h + loy - 2miy + 2moy'^ - 3niy^ + imy^) + O (A) , 

(B18) 



2(A:o + Zo2/ + TOoy2 + noy3) +0(A). (B19) 



Due to the no slip boundary conditions on the surface 
of the swimmer we have = Iq = li = 0. The 
far field boundary condition on the outer velocity field 
v|r-).oo = —Ue.z yields Z — ~U/2. Matching the low- 
est order velocities through vj^ (y — >• oo) = v° (y — > 0), 
as done for the concentrations above, we find 



= - i (2too y + 3no y^) 



(B20) 



~2Zo -Yo + Xo = ^ - 2miy + 2moy'^ - Srtiy^ + 3noy^ , 

(B21) 
(B22) 



2Zo + 2Yo + 2Xo = 2too y^ + 2no y% 



where we have used dyh{y)\y^oo = Ki (see Eq. ((27)) ). 
From Eq. (|B20p we deduce that toq = no = 0. On im- 
posing the physical constraint that dyh{y) remains finite 
for y — > cx) we conclude that TOi=ni = Oin order to 
avoid divergence of the right hand side of Eq. (jB2ip . 
Finally, Eqns. (|B21IB22p yield together the two last con- 
stants Xo = (Ki - Uo) /4 and Yq = {3Uo - Ki) /4. This 
fully determines the lowest order velocity fields. To fur- 
ther relate the speed of the swimmer Uq to the balance of 
forces, one only needs to consider the Stokeslet v° ~ l/r, 
whose long range nature reflects the presence of external 
forces. For a free swimmer the Stokeslet vanishes and 
we therefore have Yq = 0. This condition determines the 
lowest order swimming speed 



(B23) 



The first correction of 0{\^) in a series expansion of U 
(see Eq. (pS)) ) is an extension of the scheme presented 
here. 



Appendix C: Energy balance in linear, nonequilibrium 
thermodynamics 

This Appendix substantiates the definition of the 
power input in Sec. IVII by providing more details about 
the underlying assumptions and calculations. As ex- 
plained in Sec. IVI Al we assume that an external appara- 
tus is connected to the system. The apparatus consumes 
energy and keeps the system in steady state. The ensem- 
ble of system and apparatus does not exchange matter 
but only work and heat with the external world. The 
first law of thermodynamics for the ensemble reads 



P\x\ — 



5Q 



(CI) 



where Pout is the work output of the system. We set 
Pout = because no external force acts on the swimmer. 
Within the framework of classical, linear, nonequilibrium 
thcrmodynamic n^^'^^ , the heat flow 5Q balances in steady 
state the entropy production inside the system and the 
apparatus. Wc can employ a local definition of the en- 
tropy production rate to rewrite the first law 



Pi„ = <5Q= / Tie, 



'sys 



e 



app 



dV 



(C2) 



with being the entropy production rates per unit vol- 
ume. Since the external apparatus is ideal, we have 
©app = 0. For the entropy production of the system 
we employ an established formulai^, which derives from 
a local equilibrium assumption. 



e.ys =^2ry^ : (Vv) - ^J^h^ (fM + 



ki 



(C3) 
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where Uik is the rate at which the reaction with index 
k changes the concentration of species i. For the model 
system considered in Sec. IVI[ the entropy production 
rate is rewritten by employing the Stokes equation Eq. 
dl]) and the chemical reactions defined in Eqns. P4I35|) 



insertion into Eg. ljCSp yields 



+ 5{r-l)g{6) (^a - * - Mb) • 



(C4) 

In order to calculate Pin, we assume that the chemical 
potentials jla-.i^b and jls are fixed by the external appa- 
ratus at f ^ oo. This, together with Eq. (jA6p . implies 
that the concentrations at the outer boundaries of the 
system are fixed. To facilitate reference, we note here 
the steady state diffusion equation 



v5, = 0. 



V • J, = V 



The solute flux boundary conditions at f 
V = arc given by 

^rjs I r—R — 0- 



(C5) 



R where 



(C6) 



Inserting the entropy production rate Eq. (|C4p into the 
expression for the power input Eq. ()C2p and using Eq. 
(|C5[) for the solute of type b we find 



TOsy^dF 



[V 



j6 + vcfaj ^-J +^J,V/2,]d\> 

i 

Dc , , 

—g{9)ifla - ^' - Mb) dA.r.=B. 

(C7) 



The potential 5* now only occurs in surface integrals, 
^'(f) decays quickly for f oo. Furthermore, using Eq. 
(jC6p leads to a cancelation of surface integrals at f = i? 
containing VP. We are thus left with 

J TesysdV=- J [jbV/ib+jaVMa +jsV/i,]dy 



Dc 



(C8) 



giO) (Ma - Mb) dAf^ 



The Gibbs-Duhcm equation, which is consistent with our 
microscopic model, reads = CiVjli — Vp. Multipli- 
cation of this equation by the center of mass flow v and 



i 

f Dc 

+ j {i^-a - h) dAr=R 



E ^'■j^M* + (Ma - Mb)] dAf^R, 

(C9) 



where the diffusion equations and V • v = were used 
to produce the boundary integrals. Due to the bound- 
ary conditions Eq. (|C6p . the last integral in Eq. (|C9[) 
vanishes. Since the swimmer is not subjected to external 
forces, the fluid flow contains no Stokeslet. With the as- 
sumptions in Eqns. (jA61 IA7[) the pressure does also not 
depend on local concentrations. Therefore, the boundary 
work of the pressure vanishes / vp dA,-;=oo ^ and can 
consequently be dropped in Eq. (|C9p . For our model, 
we conclude that it does not make a difference whether 
we consider a pressurized system (Gibbs free energy) or a 
system with fixed volume (free energy). The power input 
becomes 



(CIO) 



Employing the chemical potentials from Appendix |X] 
with the assumption of diluteness {vaCa + VbCb 1) and 
Eq. dMl) along with Eqns. (|C5IC6p we find 



/ ^ JiMi GrdAf^rx, = ^ / X! GrdAf^ 



Jfc (/cTlog [Ca/h] + Eq - Cfc) GrdAf^oo- 



(Cll) 



This formula is employed in linearized form for Eq. (|54p . 
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